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I. INTROUUCTIO/J 


onvlrOTMcl,T!mr.,"o it"?o''3!tai„'a T 

protons and clcctronfirtho c^roy raLo1i::^‘; 

International Sun-l;arth l-'xnlorcr fTsm^A ^ n*^ 

port is the result S a sSr^n’i ^ spacecraft. This re- 

radiation damage to the thermal coatin^ n^' tl possible 

both of which will brnlaS n spinning spacecraft, 

of 21 R]., perigee of 280 km, and iLlLiation of^28 5°^^* apogee 

the model to other high-altitudo ^ * %’lications of 

propriate orbit averaging! can be obtained with the ap- 

lo the best of our knowledge, this is thp ** 

size an overall quantitative environment of i. attempt to synthe- 

altitude spacecraft, usinc data from in ••■► w-energy particles for liigh- 

ing model environments oAhe morrono^goUc (> IM'^ovr' 

pod in the radiation belts of the Partli Co o ^ particles trap- 

the spatial cells used ore small voSes ‘of ^76), 

the geomagnetic field However of the tubes of force formed by 

plasma dictates a much groLrspaJxal^st^f ncar-liarth 

particle fluxes need to brobtaSed ‘"'>•■■■“8° 

tho plasma environment consist of Various spetial regions for 

WitHn these regio\srt\-f Hr^^HrS^rj^^ai 

plasmr\\e"t:hS!iath1drLi‘:t^"Slr°^ -gnetosheath. 

been chosen as the spatial strnctnno + 1 ^ inner magnetosphere) have 

real geometries of thes“jcgl™s . t, ‘"o 

the region., must bo used; then,' the avo^agf ch5ja«ertf 
in each region must be ascertained Tn thn ^ f ^ fluxes with- 

ment, the final objcctivrrrio Sehve Ll 'i “"iron- 

energy, spin-averaged spectra for nmtrm differential 

best available onvtronSf all ‘>'^“'1" "'e 

measurements published in tlie Htor-.t higli-altitudc satellite 

of the energy or ,™gularrJsoluH^r^rH “"siJerod, regardless 

altitude, poLr-orWtlng snaccci P ?. “?"n.nont. Data from low- 
of gualitativc data Kill 1 , 1 , referenced to"°' p^°" ''ubl i cat ions 

titati,-,. remsults. Oiily fhlxe-s llll 1 L “ ‘''t '1““"- 

will he included ‘••horf t-orm t- ^ ^ lony-tcnn or reported a;; tvpieal 

glected, unles.l the;: bill;;! V™ 

average environment. x„ attempt Kill' 1,1 ,Md;'il''gl'leT‘''''’'''i‘"'' *“ 
rov.e,. of the Physics assoeiatld Kith the“ 

llie basic approacli employed in the derivitinn nf ti, i ■ 
spectra ,s divided into two step.,. thl^'l i:.!;"!;.;:;,, 
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l)clon}>ing to o;idi of tlic five rc{>ioiis arc compiled, and a composite 
spectrum representing the environment of tliat region is constructed. In 
cases wliere there are substantial differences in data from several mea- 
surements, spectra of a high estimate, as well as a low estimate, are de- 
termined. The details will be given in Section III. In the second sten 
the composite spectra of each region are weighted according to the frac 
tion ot the orbital time tlie spacecraft is located in the corresponding 
region. As the orbits process about the north ecliptic pole, tlic weight- 
ing factor of cadi region changes. The details of the second step, using 

the first 6 months of the planned ISlili-A/-B orbit, will bo described in 
bection IV. 


To date, the abundance of the publislied data and the completeness of 
various regions cannot be viewed with the same standard. 
Bocau.e of the nature of the plasma characteristics, and perhaps because 
of the tocus of interest in the scientific community, the solar wind in 

medium IS the best documented region, while the mag- 
cosheath has comparatively few publications with quantitative data in 
absolute flux units. Accordingly, the resultant spectra of one region 

may have an uncertainty factor quite different from that of the spectra 
or other regions. ^ 

AS previously noted, the particle properties are distinctively dif- 
ferent in the various spatial regions. Often, it is necessary to use 
simplifying assumptions to derive representative spectra for a particular 
region. The data on the solar wind proton bulk flow, for example re- 
quire the computation of time- and spin-averaged fluxes. Details ’of the 
matliematics will be given in the Appendixes. 
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II. 


lllJi SPATIAL STRUCTURL FOR Tllli NliAI^-liARTll PLASMA ENVIRONMLNT 


For the purpose of flux averaging along a given orbit, it is essen- 
° ^ positions of the spacecraft with respect to those 

regions of the near-Farth environment that have distinctively difLrent 
particle populations. The International Magnetospheric Studv/Satollite 
Situation Center (IMS/SSC) has developed a Smputlr program ^ 
erge a model configuration of the magnetosphere with the cphemoris of an 
Earth-orbiting satellite (Fatt. aZ. , 1976). This program ha^dLect 
application in-4;be.-pres ent study. airect 

® report presented at the 10th ESLAB Symposium, (Vette et al 
976), eight regions of the near-Earth space were defined by the simple 
boundaries that are listed in Table 1. Fairfield’s model 

netiLuse boundary locations of the bou shock and the tig- 

netopause. Some of these regions are illustrated in Figure 1. Tliis ^ 

tion^fo^thrs^H^^^^r to provide a general frame of posi- 
tions for the study of particles and fields in the IMS. ^ 

1 • I ^Gcause of the lack of sufficient low-energy particle data to estab- 
lish distinctively different statistics for a local time effect il Se 

magnetosphere, the further separation of 
these two regions into dayside and nightside regions has been eliminated 

the midlatitudc magnetotail and the^ 
nZn ^ sheet have been combined and named the plasma sheet, which is a 

regior Ara^Jcsult^or^?^“"'‘'vT Particles found in this 

region. As a result of this slight modification, only five ronions as 

defined in Table 1 have been used; the previously notL IMS/SSC computer 
program is still applicable without change. 


in. DliRIVATlON OF COMPOSITF SniCTHA 


A. The Interplanetary Mudium 

Tlie low-energy particle environment in tlie interplanetary medium 
near the Larth is dominated by the "quiet" solar wind. Since most of 
the solar wind observations found in tlic literature are described in 
terms of fluid parameters (temperature, density, and bulk flow speed), 
they must be converted to a differential energy flux; the method is given 
in Appendix A. Short-term variations in the interplanetary fluxes can be 
event^ f? ^ transient solar activity. These include flare particle 

regions, and interplanetary shock- 
a ^ review, see Lin, 1974b). Shock-associated spikes 

with a duration on the order of minutes, will not appreciably affect the ’ 

orthf considered in this 'study. The prLence 

of the barth ^d its magnetosphere can also add to the near-barth, inter- 
planetary environment through reflection and emission of particles at the 
bow shock (A8br%dge et at,, 1968; Frank, 1970b; Lin et al\ , 1974). 

and mult ^he various interplanetary fluxes varies widely 

and must be taken into account in the construction of a useful composite ^ 
spect:mm. The electron and proton components of the "quiet" solar' wind 
are distinctly different in this regard. For the electron flux the nearly 
isotropic thermal component of velocity is large compared with the bulk 
flovv component that is commonly referred to as the solar wind velocity 
For protons, the opposite situation is obtained in that the thermal com- 
ponent is small compared to the bulk flow component. Tims, the electrons 

essentially isotropic, while the protons are highly 
anisotropic. (For reviews, see Montgomery, 1972; Wotfe, 197'’ ) Tvni- 
cally the protons arrive within a cone of 20°, centered along the bulk 

thrbuJrn"”" 1967). In addition, the magnitude of 

Tn Varies within the approximate range of 250 to 700 km/s. 

0 cimpare this flux with that from other regions of space where fluxes 

?[ux"?rorr’' Ular ,,r“an 

flux IS constructed. Iliis approach will be discussed more fully in the 

0 lowing paragraplis. Details are given in Apj)endixes H and C. 

A typical low-cnergy, solar wind election spectrum, with a density 
equa' to .) cm and a temiierature equal to 1.5 x lo^ K, is shown in Figure 
.{Montgomery et al. , 1968; Montgomery, 1972). Also shown is a portion 

over spectrum {Lin at at,, 1972), whicl. extends 

over 12 decades of flux for energies uj) to 100 keV. IMI’ o observations 

anisotropic flare packet {Frenk and 
1)7.) are also shown. The remaining si)octra from (XUl 5 {Oaileie 

* '? i 1"’ subsatellite 

1 ’ Hulisatellite {Un el a/., 1975a) were all 

obtained during various flare events. These events have been observed 
for energies as low as 1 keV. .>iiveu 
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I'rora tliir> compilation of spectra, a rcjircscntativo isotropic com))os- 
ite spectrum has been conf triictcil and is also sliown in I'iyurc 2, Without 
a statistical study of tl»c flare spectra at enerj^ies above 1 keV, it is 
difficult to determine tlieir contribution to the average environment. 
However, by extending the quiet-time observation of IMl’ b {Frank and 
Gumett, 1972) to higher energies, a reasonal)le estimate of the flare 
contributions to tlie quiet-time composite of I/ln ot at. (1972) is obtained. 
The IMP b spectrum at lower energies joins smoothly with tlic Lin composite 
spectrum; and, tlicrcfore, tliey arc used as tlie basis for the isotropic 
composite spectrum below 1 keV. Also shown is a spectrum representing 
the maximum anisotropic flare fluxes observed. 

An interplanetary medium, quiet-time proton spectrum for energies 
above 38 keV is shown in Figure 3 {L-ln et at. , 1973b) . Also shown is a 
derived, interplanetary medium, average proton spectrum based on enhanced 
fluxes associated with solar active regions and flares. This spectrum is 
three orders of magnitude below the IMP 6 proton fluxes seen propagating 
upstream from the bow shock with energies above 30 keV (Lin et at. , 1974) . 
These IMP 6 protons are observed when tlie spacecraft is on an interplan- 
etary magnetic field line that intersects the bow shock. In the energy 
range from 6 to 40 keV, IMP 4 lias observed transient events, which Frank 
(1970b) has identified as magnetic storm-associated fluxes. Tlicy were 
quasi-isotropic in the ecliptic plane but favored intensities arriving 
from the solar direction by factors of two to three. A nonstorm upper 
limit from IMP 4 for energies in the range of 11 to 18 keV is also shown 
in Figure 3. Protons propagating out from the bow shock along interplan- 
etary field lines have been observed by IMP 4 (Frank, 1970b) . These 
particles, shown in Figure 3, are quasi-isotropic over a hemisphere look- 
ing toward the bow shock and have energies in the range 0.2 to 4 keV. 
However, for energies below about 4 keV and looking in the solar direc- 
tion, the solar wind dominates. 

The proton comjionent of the solar wind can be characterized as a 
highly anisotropic beam of particles with a narrow energy spread. The 
peak flux occurs at an energy close to that given by a proton moving 
with the bulk solar wind speed. Two examples of instantaneous solar wind 
proton spectra, corresponding to bulk speeds of 300 and 700 km/s, are 
sliown in I'igurc 3. Over a period of time, tlie solar wind speed ranges 
between 250 and 700 km/s. Oerived density and temperature parameters are 
found to correlate with tlie solar wind speed (for a review, sec Wotfe, 
1972). Using observations from Vela 3A and 3B obtained during the jieriod 
July 19()5 to November 1907 (Uundhauean et at,, 1970), it is possible to 
derive a time-averaged, solar wind spectrum. This spectrum is shown in 
i'igure 3, and additional details of the derivation are given in Apjiendix 
11. Although this spectrum is derived using a distribution fur solar wind 
speeds obtained near solar minimum and during a rising iiortion of the 
solar cycle, distributions obtained during other jihases of the solar cycle 
are only expected to change the average spectrum sliglitly. The highly 
anisotroiiic nature of this spectrum makes it difficult to compare with 
the otlier, more isotropic fluxes. To resolve this difficulty, an effec- 
tive spi n-averageil , isotropic s]iectrum has been generated and is shown 





in Figure 3. An area that is oriented perpendicular to tl>e ecliptic plane 
and spinning about an axis tluit is also perpendicular to the ecliptic 
plane, will receive tlie same average flux from this spectrum as it would 
from the anisotropic, timc-avcragod, solar wind spectrum. Tlie derivation 
of this spin-averaged spectrum is given in Ap])cndix C. Similarly, a spin- 
averaged, bow shock spectrum has been constructed by smoothly joining the 
two bow shock-related spectra and then dividing tlie flux by two to account 
for the quasi-isotropic nature of these fluxes over one hcmispliere. The 
resulting spin-averaged, bow shock spectrum is sliown in Figure 3. 

A composite proton spectrum for tlie interplanetary medium near the 
Farth now can be constructed by adding togetlier the spin- averaged, solar 
wind and bow shock proton spectra. Tliis spectrum, shown in Figure 3, is 
isotropic and simulates the effect of tlie observed anisotropic spectra on 
a plane area that is spinning as previously described. Tliis is approxi- 
mately the geometry for a unit area on the side of the spinning ISF1;-A/-B 
spacecraft. 


B. The Magnetosheath 

Tlie magnetosheath is the region between the bow shock and the magne- 
topause as shown in Figure 1. The bow shock is the boundary defined ob- 
servationally by the abrupt changes that occur in the characteristics of 
the plasma and the magnetic field {Montgomery et at., 1970). The magne- 
topause is the boundary formed by tlie geomagnetic field as a standoff 
limit to the solar wind plasma from the interplanetary medium. 

In the magnetosheath, the post-shock plasma from the solar wind flows 
around the magnetosphere and is usually turbulent {Hwidhausen et at,, 

1969; Saudder et at,, 1973). For the low-energy particles in tliis region, 
the presentations of information in the literature arc generally divided 
into two groups: (1) in the form of differential spectra constructed 

directly from the observed fluxes (e.g», Frank, 1970a; Rosenbauer ct at,, 
1975); and (2) in terms of particle distribution functions and moments 
expressed as conventional fluid parameters that are also derived from the 
observed fluxes (c.g,, Ihmdhausen , 1970). l^ien available, data obtained 
in absolute flux units will be used for analysis. 

Figure 4 shows the electron spectra measured by the Vela 4B, IMI’ 5, 
and Apollo 14 spacecraft with the source publication listed in the figure 
legend. Tlie Apollo 14 data are given to show tliat tlie electron sjiectra 
in the distant magnetoshe;itli (~60 lip) are similar to those observed by 
Vela 4B at about 18 li];. 

Because of the lack of com])arable statistical information about tliese 
measurements, tliore is no simple way that the various spectra can lie ]irop- 
erly weighted to construct a single, smooth spectrum, i'or tlie present 
purpose oi making a crude, quantitative projection of the environment, an 
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ui)pcr-l)ouiul Mild a lowcr-houml envelopo arc constructed as tlic composite 
spectra for a hi};h and low estimate. As sliown liy the dashed lines in 
I'ijiure -1, tlie composite s])cctra are jjenorally hroador tlian any of tlie 
individual, measured siiectra, (liven tlie uncertainty as represented hy 
the raime oi variations in the estimates, tlie comjiosite spectra can be 
reiiarded as represent iny. the lony-term, averaye fluxes in tlie reyion. 

As is the case for electrons in the inter])lanetary medium, the any.u- 
lar distribution of the maynetoslieath electron fluxes shows rcl itively 
little variation (within a factor of two) with respect to the bulk flow 
direction {Mont<jomcry et at,, 19701. Therefore, no sjiin averaying is 
needed for the present flux estimates. 

Figure 5 shows the raagnetosheath proton spectra observed by the 
llliOS 1, llbOS 2, IMP 4, and IMP 5 satellites. For the low-energy protons 
in this region, the bulk flow velocity is a major contributing factor in 
the observed kinetic energies of the particles Qhmdliausen et al. , 1969). 
The proton flux, therefore, is highly anisotropic with the peak flux in 
the direction of the bulk flow as shown by the IIHOS 2 measurements at 
three azimuthal angles in Figure 5. However, in comparison with the pro- 
tons in the interplanetary medium discussed jireviously, the protons in tlie 
magnctoslieatli generally have a lower bulk flow speed and a liigher mean 
tliermal velocity. Consequently, tlie angular distrilmtioii of tlie low-energy 
protons is less anisotrojiic in the magnetosheath than in the interplanetary 
medium. In addition, proton spectra in the magnetosheath show considerably 
more deviation from a Maxwel 1-Boltzmann distribution than in the case of 
the interplanetary medium {Fuinnaano et al. , 1973). 

A. 'reviously described for the case of electron spectra in the mag- 
netosheatn, an upper-bound and a lower-bound enveloiie arc constructed to 
serve as composite spectra reiircscnting a high estimate and a low esti- 
mate in the jiresumed flow direction. Because the com]iositc spectra arc 
generally broader than any of the individual, observed spectra, the result 
of the estimated flux, when integrated over all energies, will be conser- 
vatively higher than the average situation. However, the differential 
flux from the composite siiectrum at a given energy will probably be com- 
parable to the average measurement within the range of uncertainty as 
indicated by the liigh and low estimates. 

In order to calculate the averaged proton fluxes received by a unit 
area of a spinning surface with tlie sjiin axis perpendicular to the eclip- 
tic plane, botli the flow component as well as tlie tliermal component of 
the particle velocity must be considered in the angular integration and 
the spin-averaging, processes. In practice, as descrilied in detail in 
Apjiendix I), it is convenient to make the calculation usiii)', an analytic 
distribution function that fits the composite s]H'ctra. 

In guneral, (he K-di st r Unit i on function with K - 2 (/'.a'ri/a.m, ■ ,7 a/., 
1973) is found to give a reasonable fit to a tyiiical, observed proton 


.spectrum in tlie m.-ijiiietoslicath. In the ])ro'ient analysis, the composite 
proton spectra of the maj'netoslieatli can he well aj))m).Ni mated l>y tlie K- 
distribution witli a hulk flow speed of 2'<0 km/s, tlieniial energy (k'l') of 
400 eV, and a density of eitlier 17 0 . 111 “’. or 2 cm"’ for tlie high and low 
estimates , respectively. 

The results of the calculation (sec Ajijiendix I) for details) are 
presented in figure 6. Tlie effective isotropic s])ectra shown in tlie fig,- 
ure give the equivalent isotropic environment that the jireviously sjieci- 
fied surface area will encounter in orbit. They can be used to compare 
with other isotropic spectra for the purpose of making flux estimates on 
this area, but they cannot be used for other purposes such as calculating 
an omnidirectional flux. A procedure for correctly deriving an omnidi- 
rectional flux spectrum is included in Appendix U. 


C. The Inner Magnetosphere 

The inner magnetosphere, shown as region C in Figure 1, is composed 
of the entire dayside magnetosphere and that part of the nightside magne- 
tosphere where the absolute value of the X-coordinate is less than 10 Up 
in the geocentric solar ecliptic (GSh) coordinate system. A further re- 
striction for this study is the consideration of only near-equatorial 
latitudes. This large, asymmetric volume, whicli includes the inner and 
outer radiation belts, is known to contain a sea of protons and electrons 
with energies in the range 100 eV to 50 keV (cf. , Vasyliunas , 1972). 
Observations of the elections are considerably more complete tlian those 
of the protons, but even they are lacking in spectral, temporal, and spa- 
tial coverage. This is in part a result of background problems, some of 
which arise from the presence of higher energy particles and solar ultra- 
violet radiation [Gvinyauz, 1969), and some not yet clearly identified 
{Lyons and I'iiltiams , 1976). Nevertheless, the general jiicture for elec- 
trons is one of considerable spatial structure and temjioral variation in 
the flux of tliis region. (For reviews, see Grinyauz, 1969; Vaoytiwias y 
1972.) 

It is more difficult to characterize tlie proton population in the 
inner magnetospliere since tliere are few observations covering tlie more 
distant parts of tlie dayside region lieyond an 1. (Mcllwaln parameter) of 
about 7. 0(10 3 observations {Frank and OtJans , 1970) have provided a sur- 
vey of low-energy jirotons in the evening to midniglit quadrant. (For 
reviews, see Frwik, 1970; Vasy'Hunas , 1972.) I’rotons In the plasma sheet 
continue across the inner boundary of the plasma slieet up to the plasma- 
sphere. Tlie evening plasma sheet boundary, as identified in the electron 
flux, is not observable in the jiroton flux. These jirotoiis .ire most in- 
tense in the evening sector, forming the (|ui et-t 1 me , proton ring, current 
centered near I. = (>.5. During storms, the intensity of these protons 
greatly increases, forming the storm-time riii)', current (cf. , Ft^mi:, H)()7I>)« 
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In fidilition, the tluyside majjnetospliere boundary (iitaynetojiause) has 
been observed to contract sJy.nificantly inward duriny. L-u^.c yeoiitaynetic 
storms. The yreatest observed contraction to date was associated with 
the storm of August 2, l‘)72, when the may.netopausc was rei)ortcd to be 
inside the iixi)lorer 45 satellite position at L = 5.2 {Hoffman cl at,, 
li)75). iiie eilects of sucli a contraction on the iiarticles normally re- 
siding outside the synchronous orl)it (h - (,.(,) l,uve not yet been observed 
but are certainly significant. 


Vaoyliunas (1972) has described tlie spatial structure of tl>e electron 
fluxes with energies from 100 eV to a few kcV. Tliese electrons arc found 
within the equatorial region of the maguctotail, where they form the 
plasma sheet. These intense fluxes terminate abrujitly where they define 
the inner boundary of the plasma sheet in the evening side of the magneto- 
sphere. This occurs at about 11 Rg during "quiet" times. Tliis boundary 
approaches the plasmaspherc (~6 Rg) near the midnight meridian. Moving 
from evening to earlier, local times, this boundary remains at ~11 Rg and, 
finally, intersects the magnetopause near noon. Weak electron fluxes are 
observed between the plasma sheet and tlie plasmaspherc in the evening 
sector. They are observed to be an order of magnitude or more below those 
in the plasma sheet at a few keV {Gchield and Frank, 1970). Tliey grow 
stronger as one moves to earlier, local times and expand in spatial ex- 
tent until they fill the entire region between tlie magnetopause and tlie 
plasmapause near noon. In the morning sector, tlie fluxes show a strong 
radial gradient, increasing toward tlie magnetopause. In the predawn 
Iiours, tlie fluxes become very intense and extend all tlie way from the 
magnetopause in to the plasmapause. At tliesc local times, tlicy are an 
order of magnitude larger than at other local times. Tlicsc fluxes must 
then decrease as one moves toward mianiglit so as to merge with tlie plasma 
sheet. During substorms, the plasma sJiect boundary moves toward the 
Earth in the evening sector. It is, tlicrefore, apparent that the inner 
magnetosphere will also contain appreciable plasma slicct electrons tliat 
must be included in the generation of a composite electron spectrum for 
this region. 


figure 7 is a representative compilation of the differential elec- 
tron spectra available in tlie literature. OCIO 1 {Vcaujliunaa , I9(i8) and 
OGO 3 {Irank, 19G7c; SchieUL ami Frank, 1970) spectra were obtained in 
the evening and midnight sectors, respectively, 'liie OGO 3 spectra in the 
plasma sheet at L = 9.8 and in the edge of the pliisma sheet at 1. = 9.25 
are more than an order of magnitude higher than the electron trough spec- 
trum at L = 8.8 for energies of a few keV. At L = 5.3 (in the plasma- 
sphere) and at 1, = 3.9, the spectra are between those at 1, = 8.8 and b = 
9.8 for energies of a few kcV. At the lower energies, they continue to 
rise, while the plasma sheet spectra appe/ir to fall. Tlie OGO 1 spectrum 
was obtained during a substorm and shows tlie tyiiical jilasma sheet spectrum 
extending over the evening region from 7 to 11 Ri;. At this time, the 
5-Rii spectrum has the appearance of the OGO 3 trough spectrum. The Hawk- 
eye 1 speictrum at 1. — 5.8 {humcL !■ Kind i't‘a.nk, 197(>) and the lilectron 2 
spectrum at b •= 8 {Vcruin^ cl al. , 19g(>) were botli obtained in the moniiiig 
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observations at L = 3.9 and 6 discussed earlier. These differences are 
assumed to be indicative of the temporal variations. 'Die lixplorer 45 
observations at L = 4.25 iSmith ot at,, 1976) and L = 4.6 {Wittimn nt 
at., 1973) correspond to relatively "quiet" conditions. The perpendicu- 
larly mirroring flux at L = 4.25 was about a factor of 25 higher at 30 >eV 
during the main phase of a geomagnetic storm. Tlie two ATS 5 spectra (pe- 
Foreat, 1972) at synchronous altitude (.L 6.6) correspond to near- loss 

cone and near-equatorial ly mirroring fluxes in the evening sector. The 
similar magnitudes of these two spectra indicate a nearly isotropic flux. 
As was the case for the electrons, the error in assuming isotropy will be 
small compared to the other uncertainties. 

A single composite spectrum was derived from these observations and 
is also shown in Figure 8. This spectrum follows the upper range of the 
quiet-time observations for energies above 1 keV. Below 1 keV, the spec- 
trum follows the rising trend of the spectra at L = 3.9 (OGO 3), 4.25 
(Explorer 45), and 4.6 (ExiJlorer 45). This should provide a. conservative 
representation of this energy range, especially in view of the divergence 
of trends below a few keV. 


I). The Plasma Sheet 

In the model configuration of the near-Eartli space described in Sec- 
tion II, the plasma sheet is a nightsiue region centered about the mid- 
plane of the magnetotail with a thickness of 12 Rg. At a radial distance 
of 18 Rg frem the Earth, Vela satellite measurements have shown the same 
thickness for the plasma sheet near the dawn and dusk edges of the mag- 
netotail (^Akasofu et at., 1973). However, the Vela satellites observe 
that the apparent thickness of the region gradually reduces to about 6 Rg 
near the local midnight meridian. Since the plasma intensity in the ad- 
jacent high-latitude magnetotail is comparatively much lower than in the 
plasma sheet, an overestimate of the size of this region will give a 
conservative estimate of the fluxes in this part of space. 

Figures 9 and 10 show the low-energy proton and electron data record- 
ed by the Vela 3B, Vela 4B, and OGO 3 satellites during the 1965 to 1968 
time period. It is noted that the Vela 4B spectra in Figure 10 are com- 
posites that are based on a set of five spectra published by Hones et at. 
(1972). All the spectra used in this study are classified as typical 
plasma sheet protons or electrons by the experimenters. Effects of sub- 
storms and isolated localities of acceleration or injection arc regarded 
as transient phenomena and have been tentatively excluded. 

It can be seen in Figures 9 and 10 that tlie spectra, as observed by 
different satellites or at different times, arc quite variable in energy 
range as well as in flux level. Similar to tlie case of tlic magnetosheath, 
liigh and low estimates arc constructed to serve as composite spectra and 




11 


arc shown by the dashed lines in i'ii’ui'cs 9 and 10. Tlie area enveloped 
l)y the composite s])ectra represents tlie ranp.c of uncertainty in the esti- 
mates , 

'Fliere liave been reports that tlie protons in tlie plasma slieet form 
enhanced flows (e.g., Honco cb al., 1972, 1973, 1974) during magneto- 
spheric substorms and sometimes produce explosive jet streams (called 
"fireballs") in some localities in the magnetotail {Fvank et al,, 1976). 
However, with the exception of these occasional events during substorms, 
the fluxes of low-energy protons and electrons in the plasma sheet have 
been characterized as "most often isotropic" {Akasofu et at., 1973; Bame 
et at., 1967). Because only long-term effects are of major concern in 
the present study, no spin averaging is required, and all fluxes are taken 
to be isotropic. 


E. The High-Latitude Magnetotail 

The plasma intensity in the high-latitude magnetotail is generally 
very low during "quiet" times. It is usually below the detection thresh- 
old of the instrument or below the cosmic-ray background level. The Vela 
electrostatic analyzer, for example, can only measure the electron flux 
in the region, but cannot record the proton flux during nonstorm time, be- 
cause the proton count rate is below the count rate produced by the cos- 
mic rays in the instrument {Akasofu et at,, 19731. Comparatively, the 
Vela instruments have documented very extensive data of both electrons 
and protons for the plasma sheet region previously discussed. 

In fact, there are only two spectra published in units of differen- 
tial flux available in the literature for the high-latitude magnetotail 
region: one for protons from an IMP 4 measurement {Fvank, 1970a), and 

one for electrons from Vela 5B, 6A, and 6B observations {Akasofu et al, , 
1973). Tliese two spectra, shown in Figure 11, have already been quoted 
as representative spectra in many review articles (e.g., Vasyliunas , 1972; 
Fvank, 1970a; Wolfe and Tntviligatov , 1970). 

The shapes of the spectra in this region arc similar to those found 
in the magnctoshcath, but tlie flux level is much lower than in the mag- 
nctoshcath. The high-energy skew of the proton spectrum above 2 keV 
resembles in shape the skew of the reported ring current protons. (I’or a 
review, see Vasyliunas, 1972.) 

Isotropic angular distributions for the fluxes of protons and elec- 
trons in this region are reported as a usual feature {Fvank, 1970a; Aka- 
sofu ct al. , 1973); and, therefore, no spin averaging is needed. 

Tlic validity of these spectra in tlie liigli-latitude magnetotail is 
limited to a radial distance below 12 R|.;. 'Ihe higli-latitude magnctotail 
at lunar distances of 60 R]-;, for instance, can have quite different 
cliaractijr i st i cs {llavJy <’/, a!,, H>7(i), 
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treated separately in the p^LLt^tudf >^-cn 

layer have been observed to i-i t^* particles in this bound-irv 
magnetosheath and the plasma sheet Since'^tv between that of the 

than 1 He in most cases, it is f f! is very thin, less 
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F. 


A Su-ary of tha Composite Proton and Electron Spectra 

* 


Composite proton spectra from rw«a ^• 

12. High and low estimLes are giv-^ i” Figure 

spectra and plasma sheet spectra^ Tho^n spm-averaged, magnetosheath 
envelopes of observed differential fi ^gh-estimates, corresponding to 
upper bound to the integral fJ^Ls Sm?l’ 1'' ^ overestimatf or 

a lower bound for the integral flu:^es estimates yield 

proton spectrum has the highest fluxes -F ^gh-estimate magnetosheath 
to 6 keV; while above 6 keV the energies in the range loo eV 

The spin-averaged, interplanetary Jrotorfpecr^''®'’\®^®‘'^^^ dominates, 
time-averaged, solar wind spectrL^ shoJs « 1^™’ includes the 

ranp from about 300 eV to 3 keV *Th^ h- ^^^Fuficant fluxes in the energy 
With energies above approximately 4 keV interplanetary tail, 

upstream from the bow shock its Ltf ^ ‘^°^tesponds to protons propagating 
location and interplanetar/mag^etirfie?H the sateHUe""' 

spectra tend to be an order of magnituHr. orientation. The plasma sheet 
tosheath, while the high-latitudf.^m‘^ or more below those of the magne 
Of masnitude bel™ thofc In JL p^SrshL"! 

Ir should be noted t'h'ii- •f-i. ^ • .. 

are not true isotropic fluxesj therefore magnetosheath spectra 
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in .■iSrr^i?^ - ^'-wn 
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energies above 50 eV, while the high-i ^i’oetrum for 
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regarr'to'?is'L'd^a?S'd“a«ge°%\:!;u tegion, in 

As an example, Section V describe- th^ ’ depends on the satellite orI)it 
the Isiii:-V-B orbit. the application of this environment J,'. 


IV. 


THE; ISIili-A/-B ORBIT 


As introduced earlier. Figure I shows the projection of the model 
magnetosphere and some of the revolutions of the ISE0-A/-B orbit rotated 
into the X-Y plane of the GSE coordinate system. Tlie schematic defini- 
tion of the five spatial regions is also shown in the figure. Because 
the Z-axis is not shown, regions D and E are represented by the same area 
the X— Y projection. Of particular interest, in this case, is the evo- 
lution of the orbital position with respect to the various spatial regions. 
Starting from launch on Day 288, 1977, the spacecraft crosses the bow 
shock until late January 1978. It enters the high-latitude magnctotail 
in late Feb^^ary of the same year. About 6 months after launch, the apo- 
gees of ISEE-A/-B process 180* from the dayside to the nightside of the 
magnetosphere. 

Figure_14 shows the details of the first 6 months of orbital time 
in which IS3E-A/-B are located in various regions. It can be seen that 
the major regional contributions to the orbit come first from the inter- 
planetary medium; then from the raagnetosheath plus the inner magnetosphere; 
and, finally, from the plasma sheet and the high-latitude magnetotail. 

When the orbit-averaged spectra are calculated, the fraction of orbital 
time in each region is used as the weighting factor for that region. Re- 
sults of these computations are presented in Section V. 


As previously noted, the classification of spatial regions is limited 
by the constraints of the current computing programs and the availability 
of the low-energy particle data. Some areas of relatively less extent 
but of considerable scientific significance, such as the plasma mantle 
or boundary layer, are omitted in the present version of the environmental 
estimate. 
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Till: APPLICATION OP' THE ENVIUOlMiNT TO Till: IS1:1:~A/-B ORBIT 


The number of days per revolution in a given region for each of the 
first 76 ISE12-A/-B revolutions is shown in Figure 14. Tlus group of revo- 
lutions comprises approximately 6 months and corresponds to an apogee pre- 
cession from late morning to near midnight. Tins sample of revolutions 
can be used to generate an average revolution in which the spacecra t 
spends a certain percentage of its time in each of the five regions. 
Weighting each region spectra by the corresponding percentage of time 
spent in the region and summing the resultant spectra yields an orbit- 
averaged spectrum. Figures 15 and 16 show the results of this analysis 
for the ISEE-A/-B orbit as applied to the proton and electron spectra, re- 
spectively. The percentage of time spent in each region is indicated in 
both figures along with the weighted spectra and the final orbit-aver- 
aged spectrum. Most of the time is divided among of the inn^r 

magnetosphere [32 percent), interplanetary medium (^5.2 percent), and 
magnetosheath (20.8 percent); while little time is spent in the plasma 
sheet (12 percent) and high-latitude magnetotail (10 percent). In order 
not to underestimate the final fluxes, the high-estimate spectra were 
used when available. 

It is apparent from Figure 15 that the magnetosheath is the major 
contributor to the average proton spectrum for energies below 5 keV, while 
the inner magnetosphere is the major contributor at higher energies. If 
the low estimate had been used for the magnetosheath, the average spectrum 
would have been reduced less than a factor of three for energies below 
5 keV and an insignificant amount at higher energies. 

The corresponding analysis for the electron spectra is shown in Fig- 
ure 16. Approximately half of the average spectrum is derived from the 
inner magnetosphere. The remaining half is contributed by the high-esti- 
mate magnetoshcath spectrum for energies below 300 cV and ''JS;- 

estimate plasma sheet spectrum for energies above 300 eV. If the 1 
estimate spectra for the magnetoshcath and plasma sheet had been used, 
the average electron spectrum would liave been reduced less than a factor 

of two. 

The significance of the inner magnetosphere proton and electron 
spectra for energies in tlie range 100 eV to 100 keV is now apparent. As 
discussed in Section III, this region is characterized by large ^P^tial 
and temporal variations for low-energy ])articlcs. It is likely th.. . the 
■ompositc spectra for this region overestimate tlie average situation, but 
it is difficult to say by what factor. More observations, on the d.iy. u c 
particularly, are needed. 

The orbit-averaged proton and electron differential energy s)^ 
from Figures 15 and K> are compared in Figure 17. The diHerential flux 
|protons/(cm'’‘s*eV) I tlmmgli a unit area, oriented perpeiidiuilai to the 


ecliptic plane and spinning, may be obtained by multiplying tlie spectra 
in Figure 17 by a factor of ir. If the radiation d':,..agc expected from 
these particles is not strongly energy dependent, it is useful to have 
integral flux values. These values have been calculated by integrating 
the spectra of Figure 17 from a given energy up to 10'’ eV and then multi- 
plying by IT to convert to units of integral area flux (Ja) as shown in 
lUgure 18. An omnidirectional electron flux can be obtained by multi- 
plying Ja by a factor of four. An omnidirectional proton flux cannot be 
calculated directly from Ja, because the solar wind and magnetosheath 
proton spectra are highly anisotropic. Appendix C gives the procedure 
for an analytic calculation of the omnidirectional flux when the particle 
fluxes can be represented by a Maxwell-Boltzmann distribution moving with 
a bulk speed. A similar procedure for use with a K-distribution is given 
in Appendix D. As an aid to the user. Table 2 contains numerical values 
for the spectra shown in Figures 17 and 18. 
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Appendix A. The Conversion ITom lauid Parameter s to Hner|.y Spectra 

It is sometimes possible to represent a particle population by an 
'in'ilytic function employing only a few parameters that have ^ 

i ificancc. One such function is the Maxwell-Boltzinann distribution 

which is characterized by parameters of density (n) and temperature C) 

1953). n,. Maxwell-Bolumann phase space density distribu- 

tion function £m can be written as 


f 


n 


( 


m 

2iTkTy 


(A.l) 


where fm is the number of particles per unit (position) 

unit velocity volume, having (thermal) speed v; m is le par i , _ ^ 

Tis Ihe BolLmann constant; T is a characteristic temperature; and n is 

the number of particles per unit volume. 

another fnnetion that has been found to be useful is referred to as 
the K-distribution, which can be written as 


iK = 


nK! 


1 


(i-i^r 


(A. 2) 


given' "in the limit of K approaching infinity, fk becomes the Maxwell- 

Uoltzmann distribution fm. FormUano etal,^ ^^ll^Prclatiof 
distribution temi)crature T in terms of h using the rclatioi 


kl' = 



(A..^) 


Ihe two distributions given by (A.l) and 
angular depoiulenco. 


(A. 2) do not have an explicit 


(1 i ven 
deti'rmi ne 


a );eneral 
tlie cor res 


jiliase 
|)oiul i ng 


space density distribution f, it is 
d i fferent i a 1 -energy un i d i rect i ona 1 


useful to 
flux j. 


lixprossing these quantities 


in terms of differential elements yields 


dN 

dx dy dz dv, dVy dv. 


(A. 4) 


i _ dN 
“ dA dt dQ dli 


(A.5) 


where N is the number of 
Vx> Vy, vz are velocity 
solid angle, and 11 is an 
ume elements to give 


particles; x, y, 2 are position coordinates* 
coordinates; A is an area, t is a time, n is’a 
energy. By transforming the differential vol- 


'•V* dVy dv, = v^ dfi dv 


(A.6) 


where 
f can 


dx dy dz = dA dz = dA v^ dt 


CA.7) 


Vj^ is the component of particle 
be written as 


velocity perpendicular to area dA, 


dN 

dA dt di2 dv 


(A. 8) 


perpendicular to dA, v,^ can bo replaced by v 
Differentiating the nonrelativistic particle kinetic energy ^ ' 


to obtain 


li = 'jiiv'- 


(A. 9) 


dli = mv dv 

and substituting into (A. 8), yields 


(A. 10) 


f = 111^ . dN m . 

V dA df d: ,|g ■ 


(A.l 1 ) 
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wl)oro j is tl.L> desired flux as ^-iven in fA si Ti 

iji tA..O. Ihe iiejieraJ result of 


j 



(A. 12) 


can now be used with 
spectrum j,„: 


(A.l) 


to give the 


isotropic Maxwell -Boltzmann 


enerj’y 


T2iTk'J')^/~ 


TFf 


(A. 13) 


anisotrojl'c KemLe^of S bulJ^ 

istribution, such as the Maxwell-Bolt-m n ^ that an isotropic 

formed to a reference frame moving with tl o ion, trans- 

reasonable representation of these particle f ' velocity, provides a 
tions, the particle phase space dens it^ iistr ;''"!:'' 

forming the isotropic distribution frw^ i-s found by trans- 

with a bulk velocity Vs. Using the LiLvil?e^rr^'^^”'^^° system moving 
Phase space density remains -^tL^^inX 


where 


fi(w) = fz(v-Vs) 


(A. 14) 


w = V - Vs 


(A. 15) 


and w is tlic tliermal velocity v i< i 

IS the bulk velocity, and fj £uid f, are tl I’article velocity, V 

Thur'r stationary and moving \vference''f^V'^^’'‘'^^ ^oaisity distri 

Ihus referring to (A.l), a Maxwell^iolt-, ’ respectively, 

a bulk velocity V^ can be written distribution moving with 


t' = 


” (-A-f 

\2irkT / 


e 2k r 


(A. I(, ) 


ind using (A. 12), 


the flux is given by 


Jn, 


nm''-’ v ' _ ."il* v,r 


(A. 17) 


where 


|v-Vs|^ = v’ 




_ ■> 


V Vg cos 0 


(A. 18) 


and 0 is the angle between the observed ])article velocity and the bulk 
velocity. In a similar manner, the K-distribution flux can be found for 
a moving reference frame by combining equations (A. 2), (A. 3), and (A. 15) 
rA ^ w new phase space density distribution function; then, using 

(A. 12) to obtain the differential-energy unidirectional flux. Tlie re- 
sult gives 


= 


nK!(i=rf 


K " r(K-y (K--ii)«(kT)V^ 


1 + 


2(K-%)kT. 


V-Vc 


(A. 19) 


where |v-Vg| is as given by (A. 18). 


All calculations for this study have been made using the following 
physical constants: proton mass (mp) = 1.67 x lO'^^ kg 

electron mass (me) = 9.11 x IQ-'i kg (5.66 x 


10“^*’ eV*s»cm 
eV*K"^). 
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); Boltzmann constant (k) = 1.3S x 10"23 j.k'I (8.625 x 


20 


Appendix B. p,,„t„„ sp ectrum 

Typienl s..];ir wind fluxes, observed over ;i period of nionths. mny 
sliow hulk speed v.iri. -It ions of more tlum a factor of two. Ilsln,; a flaxweiJ- 
loltaiiiann distribution function in a reference frame iiiovin}> with the Inilk 
speed to rein-esent Vela oA and oB data, ilmulluuuu-n cl, ai, (l<)7()J have 
made a statistical study of tlie relatl onshiiis amoii}. various solar wind 
imramcters. These observations, covering the time interval July l‘.)o5 to 
.November 19(>7, were made near solar minimum and durin}> a ris inn* port ion 
of the solar cycle. Based on an 11-year cycle, the next corresponding 
period (U7() to 1<)78) sliould overla]) the LSi;ii-A/-B launch date. In any 
event, a strong solar-cycle dependence for the solar wind parameter rcla- 
tionshi))s IS not expected, 

I able 3 presents tlie Vela 3A and 3B correlations between the solar 
wind velocity (V.s) and density (n) and temperature (T) . Also given is tlie 
normalized irrequcncy of occurrence (prol, ability) of solar wind velocities 
in specified ranges. By combining tliese data with tlie exjiression for a 
.-laxwell- Boltzmann distribution in tlie moving frame as given by eipiations 
(A. 15) and (A.K,) of Appendix A, it is ])ossil)le to construct a time-aver- 
aged, solar wind proton spectrum (js). Writing the solar wind proton 

spectrum (j.s) to show the explicit deiiendcnce of n and T on solar wind 
si)cod, yields 


,i. (V,) = 


111 '^2 V' 


l2-nkT(V;jjF 


. 0 “ TklFD" ■ 2’ V, cos H) 


(B.J) 


where parameters not defined here are as given in A])pciidi 
averaged spectrum (j.s) can now be obtained from 


X A. The time- 


J» 


/ .is (V,) 


di_ 

d'Vs' 


d\' 


(B.. 


where 


/ 


df 

dV'. 


tlV\. 


(B. 3 ) 


I IS defined as the probal>ility that the Inilk solar wind speed is V.., 
and \ j,, and are tiu’ miiiimum ,unl maximum lailar wind vi' I oc i t i es pre- 

■uaited 111 lable The result in, i; speetrum, lookin,,; into tiu- bulk flow 

( 0 ), has :i imieh broader eiierg.y '.pre.id tli.in anv iiidividu;il obsnv.ition 


/ 


and is sliown in Figure 4. Tlu' peak flux occurs near 600 eV and nearly 
all the flux lies between 300 and 2000 eV. A toelmiciuc for obtaining a 
useful average spectrum from this higlily anisotropic flux is given in 
Appendix C. 
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A|)poiulix C. Thc> Deri vatioii of 


Spiii-AveraficHl, .Soliir Wiiul l^l^()t()ll Spectriiiii 


Tlu- hinl'.ly aiusotropic soJar wind ])raton flux can be reasonably an- 

distribntioi. ij, a reference frame mov- 
inj. with the Inilk Uow .sjioed. From Aiipeiidix A, t)ie differential eneri-y 
spectrum (js.) can be written as 


J, 


v' , Trr'-'-v^2v V, cowy) 

(2;ikT)yj ^ 


(C.i) 


flol? V i' particle velocity, V, is the solar wind (bulk 

flow) velocity, 0 is the anj-le between v and Vy, and kT is a character- 
istic thermal energy. Tiie otlier parameters are as defined in Apiiendix A 
lypical values for the bulk flow energy (Hf) , given by 'Mjcnuix a. 




(c:. 2 i 


1 -^ n>a«nitude larger than the characteristic tliermal energy 

(kl). Ihis anisotropic flux cannot easily l,e used for this study in the 
form given by equation (C.l). Since one objective of this study is to 
assess the radiation flux received by tlie spinning ISi;i:-A/-H spacecraft 
Jt IS appro])riate to consider the specified geometry. 

Figure 19 shows a scliematic representation of a jilane area (dA) as 
an approximation to tlie side of tlie spinning lSlii;-A/-H spacecraft, where 

b tlie ecli]itic (X-Zj plane and to tlie 

Ik flow direction (V^^). llie radiation received by this plane area will 
with tl,c anslf n a, the spacecraft spins In the presel.ee of an an U 
sotropic flux like the solar wind. Tlie ajiproadi adopted is to construct 
an effective isotropic flux (j^c) that will provide the same racnation 
I ough the plane area (dA) as the average received during one rotation. 

ether fso^r'^' 1’^' compared and combined with 

other isotropic fluxes, with the restriction that it only be used to de- 

erninie the average radiation on a spinning plane area when the spin 
axis IS perpendicular to the bulk flow. 

"'cll known, the flux through an area from one side (hereafter 
le erred to as area llux. in the presence of an isotropic difleren- 
tiaJ (lux (j) IS independent of orientation and is given by 


Ja 


/i ‘1 


✓ 


(f..i) 


wl ore dU IS a chf orential solid anj-lc, ami ty])lcal units of j. arc [j,ar 
tielcs/(cin ‘s'cViJ. It is now necessary to calculate the S])ln-averai-ed 
area flux (j^J in the solar wind. IJsiny. equation (C.l) and referring to 
tlic geometry of ligure 1‘), for a given angle n can be written as 

Ja (’l) = /j. (O') cos dll 4 , 


where 


cos 0 = cos 5 cos n + sin C sin n cos (p 


di2 = sin 5 <i5 dej) 


(C.5) 


(C.6) 


Tto^givr'^^’^^^^'^ averaging over the angle 


1 7T 

Ja = - f Ja (n) dn 


(C.7) 


where, by symmetry, the integration need only be carried out over a 
nair rotation. 


In general, combining equations (C.4) and (C.7) requires a triple 
numerical integration. However, an approximation is possible in tliis 

component of particle energy is a small fraction 
ot the buA flow component. Rewriting equation (C.4) with n = 0, gives 


it / 2 2tr 


5=0 ( J )=0 


• cos 5 sin 5 d5 d(|> (C.8) 


where 0=5 from equation (C.5). Performing tliese integrations ye 


ye i Ids 


i, foi = »(kT))^i /mvV.s 

’ T2ir7n)bmVs'' (TT 


1 e ^ 




bulTflnw''^ iw area that is iierpendicular to tlie 

bulk flow, by making tlie a])jiroximati on that this flow has no thermal 
componen* , tlie area flux at angle ri is given by 


J4 


Ja (ii) '■ Ja (0) n 


(C. 10) 


for |n| • it/ 2 and zero elsewhere. Substitut in}> jACn) into equation (C.7) 
yields 


Ja (n) " / Ja ((') ^os n dn ~ 

(1 Q Jl 


(C.ll) 


for the spin-averaged area flux. This approximate expression is accurate 
to better than 1 percent as indicated by computing tlie triple numerical 
integration implicit in equations (C.4) and (C.7) for tlie solar wind 
parameters given in Table 3. Therefore, the effective isotropic flux 
(jse)» using equations (C.3) and (C.ll), can be expressed as 


j.e 


ir 


(C.12) 


where 3^(0) is given in (C.9). Time averaging of jy^(o) has been performed 
using the procedure given in Appendix 13. iiic resulting time- and spin- 
averaged effective isotropic flux is shown in l-igure 4. 

Otlier geometries, besides that given in Figure 19, may be useful in 
some cases. One of these is the radiation received by a sphere of unit 
cross section, which requires a determination of tlie omnidirectional flux 
(jo). I'or the solar wind, jo is given by 

Jq ~ y (c 131 


where js(9) is given by equation (C.l), and dS) is given by eiiuation (C.6) 
with 0 = C* The result gives 


Jo 


CirkT 

mvVs 




(0=0) 


j. 



(C.14) 


where L/ple.il units of jo are [part i r1 es/ ( cm^ • s* eV) ] . This spectrum may 
also be time averaged using the procedure in A|ipeiidix 15. 




Appendix I). The Derivation of a Spin-Avci*u}^;ed, Mi4',nctu.sl>eatli Proton 
Spectrum IJsiii^ a K-l)istril)Ution ITmction 


The anisotropic magnetosheath proton flux can be reasonably apj)rox- 
imated by an isotropic distri1)Ution function in a reference frame moving 
with a bulk flow speed. The K-distribution function, with K = 2, 1ms 
been found to be appropriate in many cases {Fomiioano et at, , 1P73) and 
also approximates the composite spectra shown in I'igui^'o 7. I’rom Appen- 
dix A, the differential energy spectrum (jj() with K = 2 can be written as 

• - 4nm'^^ 1 

' "' (kT)« • ^ , V, ^ - 2v V, cos 6)]’ 


where v is the observed particle velocity, Vf is the bulk flow velocity, 
0 is the angle between v and Vf, and kT is a characteristic thermal en- 
ergy. Tlie other parameters are as defined in Appendix A. Unlike tlie 
situation in the solar wind, the bulk flow energy (Ef) given by 


E, = ‘imV, 2 


(b.2) 


is of the same order as the characteristic thermal energy (kT) . Never- 
theless, the anisotropy is generally sufficiently large that it is neces- 
sary to calculate an effective isotropic magnctosheatli proton flux (j^^) 
as was done for the solar wind protons. Tliis flux will provide the same 
radiation through a plane unit area as tlie average received by a spinning 
unit area with its spin axis perpendicular to the bulk flow velocity. If 
the spin axis is not perpendicular to the bulk flow velocity, this ap- 
proach will overestimate A detailed description of tlie geometry is 

given in Appendix C and is shown in I'igure 19. 

The flux through tlie unit area (area flux -j^^) for arbitrary angle 
n can be written as 


Ja (n) - Jjk? (^0 L dJi! 


(b.:5) 


where 


cos 0 “ cos cos n + sin r sin )i coj; vj> 


(D.4) 


!() 


dS2 = sin ^ dC d(|) 


(U.5) 


and jK2(0) is given by equation (D.l). The spin-averaged area flux 
is now found by averaging over the angle n to give 

1 ^ 

Ja = 7 / Ja (n) dn (d.g) 


where, by symmetry, the integration need only be carried out over a half 
rotation. Using equation (C.3) from Appendix C, the effective isotropic 
flux can now be written as 


k, = ^ (D.7) 

where the calculation of involves a triple numerical integration. 

'T'^e high- and low-estimate magnetosheath spectra of figure 7 cor- 
respond to jp(0=O) witli Vf = 280 km/s and kT = 400 eV. These spectra 
differ only in the density parameter with n = 17 cm“^ and 2 cm“^ for the 
high and low estimates, respectively. Tlie derived effective isotropic 
spectra are also shown in that figure. 

As discussed in Appendix C, flux estimates for otlier geometries may 
sometimes be useful. In particular, the omnidirectional flux (jo) from 
a K-distribution, with K = 2, may be readily found from the expression 


Jo - /Jk 2(0) dfi. 


where jpCQ) is given by equation (D.l), and is given by equation 
(D.5) with 0 = The result can be written as 


• [i * fr''"'''’'] - * 




1 ) 


(D.'J) 


where typical units of jo are [parti cles/ (cm^ 'S'eV) J . 
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Table 1. Spatial Regions of the Near-liarth linvironment 
for Uigli-Mtitudo Sateiiites 



night 

: Regions {Vctte ct at, , H176) 

Rive 

Regions (Present Study) 

1. 

interplanetary medium 

1. 

interplanetary medium 

2. 

dayside magnetosheath | 

/ 

2. 

magnetosheath 

3. 

nightside magnetosheath J 



4. 

dayside raagnetospliere | 

! - 

inner magnetosphere 

5. 

nightside magnetosphere / 

(-10Re<Xgsh<0) j 

1 


6. 

midlatitude magnetotail '' 

(Xgsm<-iore, 

1 



2Re< I ZGSM~Z neutral sheet | ^6 Re)| 

> 4. 

plasma sheet 

7. 

neutral sheet region | 

(Xgsm<-iore, ’ 




1 ZGSM-Zneutral sheet | -ZRp) j 



8. 

high- latitude magnetotail 

(xgsm<-iore. 

5. 

high-latitude magnetotail 


1 ZGSM“Zneutral sheet | ^bRp) 






ISlili-A/-B Orbit -Averaged ; 5 ])ectru 


Energy (eV) 

do 

Differential I’lux 
[electrons/ (cm^ • s • sr • eV) ] 

Jap* 

Effective Differential Elux 
[protons/ (cm®»s»sr»eV)] 

1 X 10^ 

1.6 X 10® 

1.2 X lO" 

4 X 10^ 

2.0 X 10® 

4.2 X lo" 

6 X 10^ 

1.2 X 10® 

4.7 X 10" 

1 X 10^ 

6.8 X lO" 

2.7 X lo" 

4 X 10^ 

9.0 X 10® 

1.8 X 10® 

1 X 10** 

1.3 X 10® 

4.7 X 10® 

4 X 10" 

7.0 X 10® 

1.7 X 10® 

1 X 10® 

1.1 X 10° 

8.5 X 10® 


Energy (eV) 

Jae 

Integral Area Flux** 
[electrons/ (cm®»s)] 

Jap* 

Effective Integral Area Flux** 
[protons/ (cm®*s) ] 

1 X 10® 

1.1 X 10® 

2.25 X 10® 

4 X 10® 

5.6 X 10® 

2.0 X 10® 

1 X 10® 

3.4 X 10® 

1.2 X 10® 

4 X 10® 

1.0 X 10® 

6.2 X 10® 

1 X lO" 

3.5 X 10® 

4.8 X 10® 

4 X lO" 

5.0 X 10® 

2.3 X 10® 

8 X lO" 

1.7 X 10® 

6.0 X lO" 

1 X 10® 

0 

0 


^The proton lluxcs have been sjvm averaged (see text) In such a way tliat 
tliey will give the average flux passing tlirough a spinning area wiieii 
the spin axis is i)erpendieular to tlie eclijitic plane. 
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Figure 9. Plasma Sheet Proton Spectra 
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